Chiral Properties of Quenched and Full QCD 
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We extend a technique for the chiral extrapolation of hadron masses calculated with dynamical 
fermions to those generated by quenched simulations. The method ensures the correct leading and 
next-to-leading non-analytic behaviour for either QCD or quenched QCD in the chiral limit, as well 
as the correct large quark mass behaviour. We find that the primary difference between quenched 
and dynamical baryon masses can be described by the meson loops which give rise to the different 
leading and next-to-leading non-analytic behaviour. 



Modern computing facilities, combined with innova- 
tions in improved actions for lattice QCD, mean that 
it is now possible to perform accurate quenched QCD 
(QQCD) simulations at quite low quark masses [|], ^, ||, 
^. For simulations with dynamical fermions (full QCD) 
the situation is much more difficult, but there are initial 
results at quark masses as low as 30 MeV |, |. The 
latter development has inspired studies of chiral extrap- 
olation aimed at using the full QCD data over a range of 
masses to reliably extract the physical hadron mass. 

In general, effective field theories, such as chiral pertur- 
bation theory, lead to divergent or asymptotic expansions 
1^, . While this raises doubts about the direct applica- 
tion of chiral perturbation theory to lattice data, studies 
of the mass dependence of hadron properties in QCD- 
inspired models |H, ||, 0, as well as the exactly soluble 
Euler-Heisenberg problem , suggest that one can de- 
velop surprisingly accurate extrapolation formulas, pro- 
vided one builds in the correct behaviour in both the small 
and large mass limits. For the nucleon (N) and delta (A) 
masses (and by extension all other baryons), Leinweber 
et al. |l2) have suggested an extrapolation method which 
ensures both the exact low mass limit of chiral perturba- 
tion theory (technically its leading (LNA) and next-to- 
leading non-analytic (NLNA) behaviour) and the heavy 
quark limit of heavy quark effective theory (HQET). The 
transition between the chiral and heavy quark regimes is 
characterised by a mass scale A, related to the inverse 
of the size of the pion cloud source. The rapid, non- 
analytic variation of hadron properties, characteristic of 
chiral perturbation theory, is rapidly suppressed once the 
pion Compton wavelength is smaller than this size (i.e. 
TOtt > A). 

It is straightforward to extend the method of Ref. 
to QQCD. One simply includes all the Goldstone loops 
(including both tt and rj') which give rise to the LNA 
and NLNA behaviour of quenched chiral perturbation 
theory (QxPT) [|l3[ In principle, the parameters of 
the chiral Lagrangian are dependent on the number of 
dynamical fcrmion flavours, Nj, but are independent of 



the masses of the quarks. This is a celebrated feature 
of partially quenched chiral perturbation theory ]T5| , . 
The extent of the Nj dependence is not well known, and 
a precise determination awaits dynamical fermion simu- 
lations with light dynamical quarks of varying number. 

Phenomenological investigations |l^ of the role of 
the pion cloud in hadronic charge radii indicate that 
results consistent with experiment can be obtained by 
adding fuU-QCD chiral corrections to the results of 
quenched simulations at moderate to heavy quark 
masses. This suggests that the size of the pion source 
is not changed dramatically in going from the quenched 
approximation to full QCD, motivating the use of a com- 
mon scale. A, in QQCD and QCD. We proceed to fit 
quenched and dynamical lattice data assuming negligi- 
ble Nf dependence in the chiral parameters and A. The 
extent to which these assumptions hold can only be de- 
termined via further dynamical fermion calculations in 
the light quark regime. 

By incorporating the chiral loops which give rise to the 
LNA and NLNA behaviour in QCD and QQCD we find 
a remarkable agreement between the fit parameters of 
each simulation. This supports our hypothesis that the 
behaviour of the source of the pion cloud within bary- 
onic systems behaves much the same in both QQCD and 
QCD. The differences between quenched and full QCD is 
primarily described by the differences in the associated 
chiral loops. Since the chiral corrections are expected to 
be larger for N and A than for others, this suggests that a 
similar technique may be applicable to all baryons. This 
would enable quenched simulations to play a more valu- 
able role, together with new experimental information 
from JLab and elsewhere, in deepening our understand- 
ing of baryon spectroscopy. 

With regard to the properties of the N and A we find a 
spectacular difference in QQCD. Whereas the extrapola- 
tion of the N mass is essentially linear in the quark mass, 
the A exhibits some upward curvature in the quenched 
chiral limit. As a result, the A mass in QQCD is ex- 
pected to be of the order 300-400 MeV above its mass in 



full QCD. The success of the extrapolation scheme also 
lends confidence to the interpretation of the A — mass 
splitting as receiving a contribution of order 50 MeV from 
pion loops in full QCD and up to 250 MeV in QQCD |9). 
The residual splitting in full QCD would then be natu- 
rally ascribed to some shorter range mechanism, such as 
the traditional one-gluon-exchange ||2^ . 

The method for extrapolating baryon masses proposed 
by Leinweber et al. is to fit the lattice data with the 
form: 



(1) 



where Sb is the sum of those pion loop induced self- 
energies which give rise to the LNA and NLNA behaviour 
of the mass, Mb- In the case of the N this is the sum of 
the processes N Ntt N and N — *■ Att — > TV , while 
for the A it involves A Att A and A — > Ntt A. 
In the heavy baryon limit, these four contributions {B 
B'tt — > B) can be summarised as: 



'BB' 



167r2/2 



G 



BB' 



dk 



(2) 



where w(fc) — ^JW^+rn^ and lubb' — {Mg' — Mg), and 
the constants Gbb' are standard SU(6) couplings fl^ . 

The factor u{k), which acts as an ultraviolet regulator, 
may be interpreted physically as the Fourier transform of 
the source of the pion field. Whatever choice is made, the 
form of these meson loop contributions guarantees the 
exact LNA and NLNA structure of chiral perturbation 
theory (^PT). Furthermore, such a form factor causes 
the self-energies to decrease as for TOtt >> A. One 

commonly uses a dipole, u{k) — (A^ — ^^)^/(A^ + fc^)^, 
with /i the physical pion mass. 

Quenched xPT is a low energy effective theory for 
quenched QCD ^ 0, analogous to xPT for fuU QCD 
pl| . Sea quark loops are formally removed from QCD 
by including a set of degenerate, bosonic quarks. These 
bosonic fields have the effect of cancelling the fermion 
determinant in the functional integration over the quark 
fields. This gives a Lagrangian field theory which is 
equivalent to the quenched approximation simulated on 
the lattice. The low energy effective theory is then con- 
structed using the symmetry groups of this Lagrangian. 

A study of the chiral structure of baryon masses within 
the quenched approximation has been carried out by 
Labrenz and Sharpe [p^ . The essential differences from 
full QCD are: a) in the quenched theory the chiral co- 
efficients differ from their standard values and b) new 
non-analytic structure is also introduced. The leading 
order form of the baryon mass expansion about m„ = 
is 
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FIG. 1: Quark flow diagrams of chiral 77' loop contributions 
appearing in QQCD: (a) single hairpin, (b) double hairpin. 



where the coefficients of terms non-analytic in the quark 
mass are model-independent. Throughout we use cou- 
plings as given in Ref. In addition, we have included 
an octet-decuplet mass splitting to explicitly give a value 
for cf]^ [|l^ . We also stress that the term in is absent 
in full QCD — such a term being unique to the quenched 
case. 

In fitting quenched data we wish to replicate the anal- 
ysis for full QCD while incorporating the known chiral 
structure of the quenched theory. The meson-loop, self- 
energy corrections to baryon masses can be described in 
the same form as for full QCD. The effect of quenching 
appears as a redefinition of the couplings in the loop di- 
agrams in order that they yield exactly the same LNA 
and NLNA structure as given by QxPT. For example, 
the analytic expressions for the pion-cloud corrections to 
the masses of the N and A have the same form as the 
full QCD integrals (Eq. 2h with redefined quenched cou- 
plings. We refer to Ref. [ IM for details. Assuming a weak 
Nf dependence of the chiral parameters, we describe the 
quenched self-energies using the same tree level values of 
D ^ 0.76 and F = 0.50 as in fuU QCD. 

In addition to the usual pion loop contributions, 
QQCD contains loop diagrams involving the flavour sin- 
glet 77' which also give rise to important non-analytic 
structure. Within full QCD such loops do not play a role 
in the chiral expansion because the 77' remains massive 
in the chiral limit. On the other hand, in the quenched 
approximation the 77' is also a Goldstone boson |2|] 
and the 77' propagator is exactly the same as that of the 
pion. 

As a consequence there are two new chiral loop contri- 
butions unique to the quenched theory. The first of these. 
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<Jg' ' , corresponds to a single hairpin diagram such as 
that indicated in Fig. |l|(a). This diagram is the source of 
the term proportional to 777^ (involving the couplings 7 
and 7' [0) in the chiral expansion Eq. (||). The struc- 
ture of this diagram is exactly the same as the pion loop 
contribution where the internal baryon is degenerate with 
the external state. The second of these new 77' loop di- 
agrams, arises from the double hairpin vertex as 
pictured in Fig. ^b). This contribution is particularly in- 
teresting because it involves two Goldstone boson prop- 
agators and is therefore the source of the non-analytic 
term linear in 777^. 

The total meson loop contribution to the baryon self- 
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FIG. 2: Various self-energy contributions to Ma for dipole 
mass, A — 0.8 GeV. Shown are those of quenched QCD: 

(a) a^A, (b) 5liv, (c) 5^''"'; and fuU QCD: (d) al^, (e) al^; 
and net totals: (f) Ea, (g) Sa- 
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FIG. 3: Fit (open squares) to lattice data ([^]: Quenched A, 
Dynamical A) with adjusted self-energy expressions account- 
ing for finite volume and lattice spacing artifacts. The con- 
tinuum limit of quenched (dashed lines) and dynamical (solid 
lines) are shown. The lower curves are for A'^ and upper for 
A. 



energy within the quenched approximation is given by 
the sum of these four diagrams: 
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As the resultant pion couphngs in QQCD are quite a bit 
smaller than the corresponding full QCD couplings, is 
smaller in magnitude than Tib- We display the individual 
contributions to the A mass for both quenched and full 
QCD in Fig. |^. It is notable that ct^^ and a^^^^ are 

repulsive, with o'^^''^'' vanishing, so that at light quark 
masses the total quenched chiral loop contribution to the 
A mass is repulsive, whereas it is attractive in full QCD. 

It is now straightforward to fit the quenched lattice 
data with the form: 



Mb = aB + PBfnl + i;B(m,r, A). 



(5) 



Once again the linear part describes how the mass of the 
pion-cloud source varies with quark mass. This form in- 
cludes the expected behaviour of HQET where the tt and 
r]' loop contributions are suppressed. Since, as discussed 
earlier, the meson-baryon vertices are characterised by 
the source distribution, which is argued to be similar in 
quenched and full QCD, we take all vertices to have the 
same momentum dependence — i.e. a common dipole 
mass, A. With this parameter fixed there are just two 
free parameters, a and /3, to fit the quenched data for 
each baryon. 

As described in Ref. |Q we replace the continuum inte- 
gral over the intermediate pion momentum by a discrete 
sum over the pion momenta available on the lattice, thus 
encapsulating finite lattice spacing and volume artifacts. 

The lattice data used for this analysis comes from a 
recent paper by Bernard et al. [p|. These results are 
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Full 


1.24(2) 


0.92(5) 


1.43(3) 


0.75(8) 


Quenched 


1.23(2) 


0.85(6) 


1.45(4) 


0.71(11) 



TABLE I: Best fit parameters for both full and quenched data 
sets for dipole mass, A — 0.8 GeV. All masses in are GeV. 



obtained using an improved Kogut-Susskind quark ac- 
tion, which is known to give good scaling properties [ p4[ . 
Unlike the standard Wilson fermion action, masses deter- 
mined at finite lattice spacing are excellent estimates of 
the continuum limit results. The physical scale of both 
full and quenched data sets has been set via a variant 
of the Sommer scale This procedure, based on the 
static-quark potential where chiral corrections are negli- 



gible 1 19 , provides a self-consistent determination of the 
scale for both simulations. 

Fits of the the form for full QCD, Eq. |l|, and quenched 
QCD, Eq. I are shown in Fig. |[ In fitting to data we 
choose a dipole mass of A = 0.8 GeV, which has been 
optimised to highlight the remarkably similar behaviour 
of the pion-cloud source in both quenched and dynamical 
simulations. Phenomenologically, this agrees with quite 
general expectations that it should be somewhat smaller 
than that for the axial form factor H^, |2^, |2^. The 
parameters obtained from our fits are shown in Table ^. 

We highlight the fact that the best fit parameters, a 
and (3, obtained from both the quenched and full simula- 
tions, agree within errors. This suggests that the quark 
mass behaviour of the pion-cloud source is quite similar in 
full and quenched QCD, lending support to the hypothe- 
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sis made in this analysis. This leads one to conclude that 
the dominant effects of quenching can be attributed to 
the first order meson loop corrections. 

We have investigated the quark mass dependence of 
the N and A masses within the quenched approxima- 
tion. The leading chiral behaviour of hadron masses in 
quenched QCD is known to differ from the full theory. 
This knowledge has been used to guide the construction 
of a functional form which both reproduces this correct 
chiral structure, is consistent with current lattice simu- 
lations and encompasses the HQET properties. The suc- 
cess of this method in the quenched case further verifies 
the importance of including meson-induced self-energies 
when extrapolating lattice results. 

We find that, although the quenched approximation 
gives rise to more singular behaviour in the chiral limit, 
these contributions are quickly suppressed with increas- 
ing quark mass. In the nucleon, the effects of quench- 
ing reduce the amount of curvature expected as lighter 
quark masses are simulated. In contrast, for the A we 
find some upward curvature of the mass in QQCD as the 
quark mass approaches zero. In addition, the A-A^ mass 
splitting increases to around 400 MeV at the physical 
point. As a consequence of this behaviour, the A mass 
in the quenched approximation is expected to differ from 
the physical mass by approximately 25%. 

Our calculations suggest that the one-loop meson 
graphs which generate the leading and next-to-leading 
non-analytic behaviour are the primary difference be- 
tween baryon masses in quenched and full QCD. Thus, 
rather than quenched lattice QCD being regarded as an 
uncontrolled approximation, we are able to make a quan- 
titative estimate of errors over a range of quark mass. It 
is vital to investigate the assumptions made with further 
dynamical fermion simulations at low quark masses to 
test the extent to which the presented results hold. Nev- 
ertheless, this discovery represents a remarkable step for- 
ward in relating lattice QCD to observed hadronic prop- 
erties. 
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